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ABSTRACT 

Nonperturbative QCD is studied with the dual Ginzburg-Landau theory, where 
color confinement is realized through the dual Higgs mechanism by QCD-monopole 
condensation. We obtain a general analytic formula for the string tension. A 
compact formula is derived for the screened inter-quark potential in the presence of 
light dynamical quarks. The QCD phase transition at finite temperature is studied 
using the effective potential formalism. The string tension and the QCD-monopole 
mass are largely reduced near the critical temperature, T c . The surface tension is 
estimated from the effective potential at T c . We propose also a new scenario of 
the quark-gluon-plasma creation through the color-electric flux-tube annihilation. 
Finally, we discuss a close relation between instantons and QCD-monopoles. 
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1. Dual Higgs Mechanism for Color Confinement 



Color confinement and dynamical chiral-symmetry breaking are quite outstand- 
ing features in the nonperturbative QCD. 1 ' 2 In particular, color confinement is 
extremely unique, and is characterized by the formation of the color-electric flux 
tube 2 with the string tension about lGeV/fm. To understand the confinement 
mechanism, much attention has been paid for the analogy between the supercon- 
ductor and the QCD vacuum. 3 Similar to the superconductivity, 4 the color-electric 
flux seems to be excluded in the QCD vacuum, which leads the formation of the 
squeezed color-flux tube between color sources. 

In this analogy, color confinement is brought by the dual Meissner effect origi- 
nated from color-magnetic monopole condensation, which corresponds to Cooper- 
pair condensation in the superconductivity. As for the appearance of color-magnetic| 
monopoles in QCD, 't Hooft 5 proposed an interesting idea of the abelian gauge fix- 
ing, which is defined by the diagonalization of a gauge- dependent variable. In 
this gauge, QCD is reduced into an abelian gauge theory with QCD-monopoles, 
which appear from the hedgehog-like configuration corresponding to the nontrivial 
homotopy class on the nonabelian manifold, tt2(S\J(N c )/\](1) Nc ^ 1 ) = Z^~ x . 

We briefly compare the dual Higgs mechanism in the nonperturbative QCD 
vacuum with the ordinary Higgs mechanism. Like the Cooper pair in the super- 
conductivity or the Higgs field in the standard theory, the charged-matter field 
to be condensed is the essential degrees of freedom for the Higgs mechanism. On 
the other hand, there is only the gauge field in the pure gauge QCD, and hence 
it seems difficult to find any similarity with the Higgs mechanism. In the abelian 
gauge, however, only the diagonal gluon remains as the gauge field, and the off- 
diagonal gluon behaves as the charged-matter field, which leads QCD-monopoles 
as the relevant degrees of freedom for color confinement. Condensation of QCD- 
monopoles leads to mass generation of the dual gauge field through the dual Higgs 
mechanism, 6 ' 7 which is the dual version of the Higgs mechanism. Thus, the QCD 
vacuum can be regarded as the dual superconductor after the abelian gauge fixing. 
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In this framework, the nonperturbative QCD is mainly described by the abelian 
gauge theory with QCD-monopoles, which is called as the abelian dominance. 
Many recent studies based on the lattice gauge theory have supported QCD- 
monopole condensation and the abelian dominance in the maximal abelian gauge. 8 
The dual Higgs scheme predicts the existence of the dual gauge field and the QCD- 
monopole as the relevant degrees of freedom related to color confinement. 3 ' 6,7 It 
can be proved that these particles are color-singlet, so that they can be observed as 
physical states. The dual gauge field and the QCD-monopole appear as a massive 
axial-vector glueball and a massive scalar glueball, 6 ' 7,9 corresponding to the weak 
vector boson and the Higgs particle in the electro-weak standard theory. 

2. Dual Ginzburg-Landau Theory and Inter-Quark Potential 

We study the nonperturbative QCD using the dual Ginzburg-Landau (DGL) 
theory, 6 ' 7 ' 10 which is an infrared effective theory of QCD based on the dual Higgs 
mechanism by QCD-monopole condensation. 6 ' 9 The DGL Lagrangian is described 
by the diagonal gluon = (A3, Ag), the dual gauge field B^ = (B^ , B% ) and the 
QCD-monopole field Xa(a = 1,2,3), 6 ' 7 

£dgl = ~[n-(3A A)]> •* (d A B)\ v - ^([n ■ (d A A)} 2 + [n ■ (d A B)} 2 ) 

3 

+ q(i jd- eH- fi-m)q + " 9?a ■ B^Xaf 2 ~ H\Xa\ 2 - v 2 ) 2 } 

a=l 

(2.1) 

in the Zwanziger form, 11 where the duality of the gauge theory becomes manifest. 
Here, n M corresponds to the direction of the Dirac string, e is the gauge coupling 
constant, g is the unit magnetic charge obeying the Dirac condition eg = 4ir, and 
e a denotes the relative magnetic charge of the QCD-monopole field Xa- 6 ' 7 The 
magnetic charge ge a is pseudoscalar because of the extended Maxwell equation, 
V • H = p m . Hence, the dual gauge field B^ is axial- vector. In the absence of 

— * — * 

matter fields, one finds an exact dual relation between A^ and B^ in the field 
equation, d A B =* (d A A). 
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In the DGL theory, the self-interaction of the QCD-monopole field Xa is in- 
troduced to realize QCD-monopole condensation. When QCD-monopoles are con- 
densed, the dual Higgs mechanism occurs, and the dual gauge field becomes 
massive, tub = V%gv. The color-electric field is then excluded in the QCD vacuum 
through the dual Meissner effect, and is squeezed between color sources to form 
the hadron flux tube. The QCD-monopole also becomes massive as m x = 2\/\v. 

As for the symmetry of the DGL theory, there is the dual gauge symmetry 
[U(l)3xU(l)8]m corresponding to the local phase invariance of the QCD-monopole 
field Xa 6 ' 7 as well as the residual gauge symmetry [U(l)3xU(l)8] c embedded 
in SU(3) C . The dual gauge symmetry leads to the conservation of the color- 
magnetic flux. In the QCD-monopole condensed vacuum, the dual gauge symmetry 
[U(l)3xU(l)8] m is spontaneously broken, and therefore the color-magnetic flux is 
not conserved. On the other hand, the residual gauge symmetry [U(l)3xU(l)8] e 
is never broken in this process. 6 

We investigate the inter-quark potential in the quenched level using the DGL 
theory. 6 By integrating over and in the partition functional, the current- 
current correlation 6,7 is obtained as Cj = —^j^D^j^ with the nonperturbative 
gluon propagator, 

if VuVv 1 1 m r 1 a ft r) 

D ^ = [9»» + («c " ~ ^2_ m | (n . p) 2 e mP^n nVp 

(2.2) 

in the Lorentz gauge. Putting a static quark with color charge 1,6 Q at x = r and 

— * 

a static antiquark with color charge — Q at x = 0, the quark current is written as 
jfi(x) — Qgfio{$ 3 ( x — r) — 5 3 (x — 0)}. Because of the axial symmetry of the system 
and the energy minimum condition. 6 , one should take n//r, which is also used 
in a similar context of the dual string theory. 3 Then, one obtains the inter-quark 
potential including the Yukawa and the linear parts, 6 ' 7 V{r) = — ^ • e r . B + kr. 

To derive the expression for the string tension k, we consider an idealized 
long flux-tube system, where the field variables can be described by the cylindrical 
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coordinate ry = (x 2 + y 2 ) 1 / 2 . Like the Abrikosov vortex in the superconductivity, 4 , 
one should consider the "core" of the hadron flux tubes, where the QCD-monopole 
condensate x almost vanishes. The analysis of the field equation shows x( r r) — 
m x vrT inside the core [tt ~ m x 1 )i an d x( r r) — v outside the core (ry ^ m x 1 ^- 
Hence, we adopt the Lorentzian-type ansatz for the QCD-monopole field, 

2 2 2 

X 2 (4)~v 2 2 or X \pt)^v 2 ™* (2.3) 

1 + m x r^ + m£ 

with p T = (p 2 +Py) 1 ^ 2 — r T l . In this case, we obtain an analytical expression for 
the string tension, 

fc= ^g^ ln X V * = 

^jm\-Am 2 B \m x - ^m 2 x -Am 2 B ) 

For the type-II limit (rng <C m x ), one finds k ~ ^g^ 5 ln(^f), 6 corresponding to 
the well-known formula for the energy per unit length of the Abrikosov vortex in 
the type-II superconductor. 4 

As for the parameter set, the recent lattice QCD studies 12 suggest mg ~ m x , 
which means the QCD vacuum corresponds to the dual-superconductor of the 
type between type I and type II. We show in Fig.l the inter-quark potential with 
the choice of e = 2.0, m x = tub = 1.67GeV corresponding to A = 29.4 and 
v = 0.154GeV, which provide k ~ 0.9GeV/fm for the string tension and the radius 
of the hadron flux as R ~ 0.12fm. Here, we have included the correction coming 
from off-diagonal gluons in the short-range part. 



' ; " ' ' ' • , , , 

_arccos — — . (2.4) 
4m 2 -ml 2mB 



^m B m x m x 
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3. Infrared Screening Effect by Dynamical Light Quarks 

We study the infrared screening effect on the confinement potential due to light 
quarks. 6 For instance, a long hadron string can be cut through the light q-q pair 
creation, and therefore the inter-quark potential is saturated in the infrared region. 
Such a tendency is observed in the lattice QCD with dynamical quarks. 2,13 

We estimate the q-q pair creation rate w in the color-electric field inside the 
hadron flux tube, which is formed between valence quarks. Using the Schwinger 
formula 6,14 in QCD, we estimate the expectation value of the energy of the created 
q-q pair as (E qq ) ~ 0.85GeV. Since the energy (E qq ) is supplied by the missing 
length of the hadronic string, the infrared screening length R sc can be estimated 
from kR sc ~ (E qq ). Thus, one obtains i? sc ~ lfm, which corresponds to a typical 
value of the hadron size. 6 

The hadronic string becomes unstable against the q-q pair creation when the 
distance between the valence quarks becomes larger than i? sc . This means the 
vanishing of the strong correlation between the valence quarks in the infrared 
region, so that the corresponding infrared cutoff, a ~ R^, 1 ~ 200MeV, should 
appear in the system. 6 Taking account of the infrared screening effect, we introduce 
the infrared cutoff a to the nonperturbative gluon propagator (2.2) by replacing 
(npp ~~ * (np^+a 2 ' 6 k ecause the non-local factor (x\ jz^p. \y) provides the strong 
and long-range correlation as the origin of the confinement potential. 6 Here, this 
gluon propagator keeps the residual gauge symmetry. Such a disappearance of the 
infrared double pole in the gluon propagator in the DGL theory can be qualitatively 
shown by considering the polarization diagram of quarks. 

Using the above gluon propagator, we obtain a compact formula 6 for the quark 
potential including the infrared screening effect by the q-q pair creation, 

f)2 p -m B r i _ p ~ar 

V sc (r) = ■ - + k ■ (3.1) 

An r a 

which exhibits the saturation for the longer distance than a -1 ~ lfm. This formula 
for the screened quark potential has been used not only for the lattice QCD results 
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with light dynamical quarks 2 , but also for the phenomenological analysis of the 
hadron decay. 15 

4. QCD Phase Transition at Finite Temperature 

We study the QCD vacuum at finite temperature 16 using the DGL theory 17,18 
at the quenched level, where the quark degrees of freedom are frozen. In this case, 
the quark term can be dropped in the DGL Lagrangian, and therefore one obtains 
a simple partition functional 17 ' 18 after the integration over the gauge field A^, 

Z[J\ = J V Xa VB, exp (i J d 4 x{C^l ch -JJ2 \Xa\ 2 } ) , 
^DGL Ch = ~\(d,B v - d v B,f + - gea ■ B,) Xa \ 2 - \(\ Xa \ 2 - v 2 ) 2 ]. 



(4.1) 

Here, we have introduced the quadratic source term. 17 ' 18 The thermodynamical 
potential is then obtained as 



V eS ( X ; T) = 3X(x 2 - v 2 ) 2 + 3 J / dkk 2 In (l - e v 7 ^!^ j 



o 

oo 



(4.2) 



+ 3 -^J dkk 2 \n (l- e -v^+^/r). 



The glueball masses m x and mg depend on the QCD-monopole condensate \i 

m 2 x (x) = 2\(3x 2 -v 2 ) + J(x)=4\ X 2 , m 2 B ( X ) = 3g 2 x 2 , (4.3) 

which are always non-negative for the whole region of X - Owing to the introduction 
of the quadratic source term in Eq.(4.1), we can formulate the effective action for 
the whole region of the order parameter without any difficulty of the imaginary 
scalar-mass problem 16 47,18 j n ^he 4 -type theory. 
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Like the Ginzburg-Landau theory in the superconductivity, one should con- 
sider the possibility of the T-dependence on the parameters (A,i>) in the DGL 
theory. In particular, the self-interaction of Xa is introduced phenomenologically, 
and it should be reduced at high T according to the asymptotic freedom behav- 
ior of QCD. Hence, we adopt a simple ansatz for the T-dependence on A, 17 ' 18 
A(T) = A(l - aT/T c ). (We take A(T) = for T > T c /a.) We take a ~ 0.97 to 
reproduce the thermodynamical critical temperature T c ~ 0.2GeV, which means 
a large reduction of the self-interaction among QCD monopoles near T c . Our re- 
sults to be shown below do not depend largely on the value of T c , which may 
takes a larger value, e.g. T c ~ 0.26GeV, suggested from the recent lattice QCD 
simulations 19 . 

We find a first-order phase transition at T c =0.2GeV. The lower and upper 
critical temperatures are T\ ow = 0.163GeV and T up = 0.205GeV, respectively. 
We show in Fig. 2 the glueball masses wb(T) and m x {T) at finite temperature. 
A large reduction of the glueball mass is suggested near T c . In particular, the 
QCD-monopole mass m x (T) largely drop down to m x ~ T c (~ 0.2GeV) near T c . 
The deconfinement phase transition occurs at the temperature satisfying m x ~ 
T, which seems quite natural because only low-lying modes with u n ^ T can 
contribute significantly due to the thermodynamical factor l/(e tJn / T ±1). Similar 
glueball-mass reduction is also suggested by the thermodynamical studies based 
on the lattice QCD data. 20 

We show in Fig. 3 the string tension k(T) at finite temperature, calculated by 
using Eq.(2.4). The string tension k(T) decreases rapidly with temperature, and 
drops down to zero around T c = 0.2 GeV. Hence, one expects a rapid change 
of the masses and the sizes of the quarkonia according to the large reduction of 
k(T) near T c . Our result agrees with the lattice QCD data in the pure gauge 21 : 
k(T) ~ k(0)(l-T/T c ) 0A2 . 

We estimate the surface tension a between the confinement and deconfinement 
phases using the effective potential at T c in the DGL theory. There are two minima 
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at x — and \ = Xc in ^eff(x> T c ). The mixed phase includes both the confinement 
phase (x = Xc) and the deconfinement phase (x = 0). When the boundary surface 
in the mixed phase is taken on the xy-plane (z — 0), the system depends only on 
the z-coordinate, and the boundary condition is given as x(z — —oo) = 0, x(z = 
oo) = Xc- The surface tension a in the DGL theory is estimated as 



oo , 2 



r ^ I ' {'[^) -M;, tt .[v(:):7;.]^. (1.1) 



We approximate the figure of V e s(x; T c ) (0 < x < Xc) as a sme curve, 
V e s(x',T c ) ~ ^{1 — cos(2nx/Xc)}, where h and Xc corresponds to "height" and 
"width" of the "potential barrier" between the two stable states at T c . The field 
equation of x(z) is then solved analytically like the sine-Gordon equation, 22 

X(,)^tan-V/*, 5=^Xc/Vh, ^—Vhxc, (4.5) 

6 71 71 

where 5 denotes the thickness of the boundary between the two phases. In terms 
of the effective potential, the smallness of a corresponds to the weakness of the 
first-order phase transition, because a takes smaller value for smaller "height" h 
or "width" Xc of the "potential barrier" in V e ff(x;T c ). 

One finds Xc — 0.75fm _1 and h ~ 0.33fm -4 from V e s(x; T c ). Hence, the surface 
tension is estimated as a 1 / 3 ~ 196MeV, and the thickness of the border between 
the two phases is 5 ~ 0.7fm. Since the above estimation has been done in the 
quenched level, the obtained results are to be compared with the lattice QCD data 
in the quenched level, e.g. o -1 / 3 ~ 80MeV. 23 Therefore, our estimation for c 1 / 3 
seems rather good in spite of the rough treatment. 
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5. Application to Quark-Gluon-Plasma Physics 



We apply the DGL theory to the quark-gluon-plasma (QGP) physics in ul- 
trarelativistic heavy-ion collisions. In a standard picture of the QGP formation, 
many color-electric flux tubes are formed between heavy ions immediately after the 
collision. 14 ' 24 In this pre-equilibrium stage, there occurs q-q pair creation violently 
inside tubes. 14,24 , and the energy of the color-electric field turns into that of the 
stochastic kinetic motion of quarks (and gluons). The energy deposition and the 
thermalization thus occur. 

We here examine the interaction between the color-electric flux tubes in the 
DGL theory to study the QGP formation in terms of the flux-tube dynamics, 
because many flux tubes would overlap in the central region between heavy ions 
just after collisions. There are several kinds of flux tubes in the QCD system. Each 
flux tube is characterized by the color charge 1,6 Q at its end. 

We study the interaction between two color-electric flux tubes with the color- 

— * — * 

electric charge Q\ and Q2 at their ends. The system is idealized as two sufficiently 
long flux tubes, where the separation distance between them is denoted by d. For 
d ^> rn^ 1 , the interaction energy per unit length in this system is estimated as 17 
E- mt ~ Q)^ 2 m 2 B KQ{rriBd) using the similar calculation for the Abrikosov vortex in 
the superconductor. 4 

There are two interesting cases on the interaction between two color-electric 
flux tubes. 

(a) For the same flux tubes with opposite flux direction (e.g. R-R and R-R), one 

— * — * — * — * cy 

finds Qi = —Q2 i.e. Qi ■ Q2 — —e /Z, so that these flux tubes are attracted 
each other, and would be annihilated into dynamical gluons. 

(b) For the different flux tubes satisfying Qi ■ Q2 < (e.g. R-R and B-B), one 
finds Qi • Q2 — — e 2 /6, so that these flux tubes are attractive, and would be 
unified into a single flux tube (similar to G-G flux tube). 
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Based on the above calculation, we propose a new scenario of the QGP for- 
mation via the annihilation of the color-electric flux tubes. 17 When the flux tubes 
are sufficiently dense in the central region after the collisions, many flux tubes are 
annihilated or unified. During their annihilation process, lots of dynamical gluons 
(and quarks) would be created, and the energy of the flux tubes turns into that of 
the randam kinetic motion of gluons (and quarks). The thermalization is achieved 
through the stochastic gluon collisions, and finally a hot QGP would be created. 
Here, the gluon self-interaction in QCD plays an essential role to the thermalization 
process, which is quite different from the photon system in QED. 

In more realistic case, both the quark-pair creation and the flux-tube annihila- 
tion would take place at the same time. For instance, the flux tube breaking 14 ' 24 
would occur before the flux tube annihilation for the dilute flux tube system. On 
the other hand, in case of extremely high energy collisions, these would be lots of 
flux tubes overlapping in the central region between heavy ions, and therefore the 
flux tube annihilation should play the dominant role in the QGP formation. In 
any case, the DGL theory would provide a calculable method for the dynamics of 
the color-electric flux tubes in the QGP formation. 



6. Relation between Instanton and QCD-monopole Trajectory 

Finally, we study the relation between the QCD-monopole and the instanton 22 , 
which is another important topological object in nonabelian gauge theories. There 
is an ambiguity on the gauge-dependent variable X(x) to be diagonalized in the 
abelian gauge fixing, and therefore we choose a suitable X(x) to describe the 
instanton configuration. The Polyakov gauge, where A±(x) is to be diagonalized, 
is very interesting, because A±(x) takes the hedgehog-like configuration near the 
well-localized instanton, and the QCD-monopole trajectory should pass through 
its center inevitably. Here, we show this relation in the Euclidean SU(2)-gauge 
theory. 
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The gauge configuration near the well-localized instanton is given by 

M*) - -^° | ( X)|° + a g (61) 

in the non-singular gauge. 22 Here, 7] afIU is the 't Hooft symbol; Xq = (xo, to) and 
ao denote the center coordinate and the size of the instanton, respectively. In 
particular, one finds 

, / x a a (yL — xn) fl 

M{x)~% K °> 2 , (6.2) 
\(x - x )\ 2 + afi 

so that there inevitably appears a QCD-monopole trajectory with temporal direc- 
tion penetrating the center of the instanton in the Polyakov gauge. For instance, 
the QCD-monopole trajectory x^ = (x, t) is exactly found as 

x = x , — oo < t < oo (6.3) 

for the one-instanton solution at the classical level. Thus, the QCD-monopole 
trajectory is expected to have a close relation to the instanton configuration. 25 

We find an interesting feature of the QCD-monopole trajectory in the Polyakov 
gauge in the multi-instanton solution 22 , 

m*) = -^* a [ e ) 1 ( x + £ ■ (6 - 4) 

For instance, there appear two junctions and a loop in the QCD-monopole tra- 
jectory in the two-instanton system as shown in Fig.4. Here, we consider the two 
instantons with the same size locating at (±xo, 0,0,0) for simplicity. In this case, 
the QCD-monopole trajectory is found to be (x, 0, 0, t) with 

x = or t 2 = (xl -x 2 ) +2\xq\sJ(xI - x 2 ). (6.5) 

The QCD-monopole trajectories tend to be highly folded by connecting their loops 
in the multi-instanton configuration. Hence, the presence of instantons is expected 
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to promote QCD-monopole condensation, which is characterized by a folded long 
monopole-loop 8 . This conjecture can be checked by the lattice QCD. 

We acknowledge O. Miyamura and H. Monden for discussions and comments. 
One of the authors (H.S.) is supported by the Special Researchers' Basic Science 
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FIGURE CAPTIONS 

1) The inter-quark potential V(r) in the dual Ginzburg-Landau theory. The 
dashed curve denotes the Cornell potential. 

2) The glueball masses mg(T) and m x (T) at finite temperature. A large 
glueball-mass reduction is found near T c . The phase transition occurs at 
the temperature satisfying m x ~ T, which is denoted by the dotted line. 

3) The string tension k(T) at finite temperature T. The lattice QCD result in 
the pure gauge in Ref.[21] is shown by the dashed curve. 

4) in the two-instanton system in the Polyakov gauge. The two instantons with 
the same size are located at (±xo, 0,0,0) shown by small circles. There 
appear two junctions and a loop in the QCD-monopole trajectory. 
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